Limits on the cosmological abundance of supermassive compact objects from a search 

for multiple imaging in compact radio sources 



P.N. Wilkinson, 1 D.R. Henstock, 1 I.W.A. Browne, 1 A.G. Polatidis, 1 ' 2 P. Augusto, 1 - 3 
A.C.S. Readhead, 4 T.J. Pearson, 4 W. Xu, 4 G.B. Taylor, 4 - 5 and R.C. Vermeulen 4 - 6 

'University of Manchester, Jodrell Bank Observatory, Macclesfield, Cheshire SK11 9DL, UK 

2 0nsala Space Observatory, Chalmers Institute of Technology, S-^3992 Onsala, Sweden 
3 Universidade da Madeira, Dep. Matemdtica, Caminho da Penteada, 9050 Funchal, Portugal 
4 California Institute of Technology, Pasadena, California, 91125 
5 National Radio Astronomy Observatory, Socorro New Mexico, 87801 
6 Netherlands Foundation for Radioastronomy, P.O. Box 2, 
7990 AA Dwingeloo, The Netherlands Dwingeloo, The Netherlands 
(Dated: February 1, 2008) 

Using Very Long Baseline Interferometry we have searched a sample of 300 compact radio sources 
for examples of multiple imaging produced by gravitational lensing; no multiple images were found 
with separations in the angular range 1.5-50 milliarcsec. This null result allows us to place a limit on 
the cosmological abundance of intergalactic supermassive compact objects in the mass range ~ 10 6 
to ~ 1O 8 M0; such objects cannot make up more than ~ 1% of the closure density (95% confidence). 
A uniformly distributed population of supermassive black holes forming soon after the Big Bang do 
not, therefore, contribute significantly to the dark matter content of the Universe. 

PACS numbers: 95.35+d 97.60.Lf 98.80.Es 98.90.+S 



I. INTRODUCTION 

The possibility that a first generation of objects with 
masses comparable with those of globular clusters formed 
prior to galaxies has long been recognised y], g. Such 
Jeans-mass (~ 10 65 M Q ) objects forming shortly after 
the decoupling of matter and radiation in the early uni- 
verse could have evolved to black holes and it is possible 
that some of the dark matter could be in this, difficult to 
detect, form || Bj. Building on these ideas Gnedin & Os- 
triker || and Gnedin, Ostriker & Rees fl explored a cos- 
mogonic model in which the baryonic density (fij ~ 0.15) 
is up to an order of magnitude higher than that in- 
ferred from primordial nucleosynthesis with the "excess" 
baryons being lost in the collapse of such Jeans-mass 
objects. It has also been conjectured that relic massive 
black holes might provide the seeds for quasars 0. 

Press & Gunn || developed the idea of detecting su- 
permassive compact objects (CO) by their gravitational 
lensing effects well before the actual discovery of gravi- 
tational lenses in 1979. They showed that in a universe 
filled with a mass density f^co ~ T the probability of a 
distant source being detectably lensed by a supermassive 
CO is of order unity, while for f2co < 1, the probability 
decreases in direct proportion to the mass density. From 
this they drew the important conclusion that the fraction 
of distant galaxies that is lensed by CO directly measures 
i7co an d, ignoring angular resolution effects, is indepen- 
dent of the mass Mco of the lenses. The latter property 
is simply understood. A given value of f^co can be made 
up of a large number of low-mass objects or a small num- 
ber of high-mass ones; hence the number density n of CO 
of a particular mass is proportional to I /Mco- For point 
masses the gravitational lensing cross-section a cx Mco 
Pj and hence the path length to lensing (1/na) is inde- 



pendent of the lens mass. However, the average image 
separation measures Mco directly and is approximately 
independent of flco- These ideas were further developed 
by Nemiroff |o) and by Nemiroff & Bistolas Q. 

Since for a lens at a cosmologically significant distance 
the image separation is ~ 2 x 10~ 6 (Mco/Mq) 1 / 2 arcsec- 
onds, searching for Jeans-mass CO requires the milliarc- 
sec (mas) resolution of Very Long Baseline Interferom- 
etry (VLBI). Kassiola, Kovner & Blandford here- 
after KKB), used the lack of "millilensed" images in pub- 
lished VLBI maps of 48 compact radio sources, to place 
an upper limit f^co < 0.4 (99.7% confidence) in the mass 
range 10 7 -10 9 M . In a closely-related search Nemiroff 
]l3[ used the lack of secondary "echoes" of gamma ray 
bursts, which would also arise from millilcnsing, to rule 
out a closure density of supermassive CO. In this pa- 
per we describe a millilens search based on high-quality 
VLBI maps of 300 sources which extends KKB's lower 
mass range to ~ 10 6 M and enables us to push their 
limit on f^co for Jeans-mass objects down by a factor 
15. This is the most stringent limit to date on flco fo r 
uniformly-distributed Jeans-mass objects and we discuss 
some cosmological implications of this result. 



II. SELECTION OF LENS CANDIDATES 

The parent sample was drawn from catalogues of com- 
pact radio sources which have been systematically ob- 
served by the authors with intercontinental VLBI arrays 
at a frequency of 5 GHz and a resolution ~ 1 mas; a 
significant fraction of these sources was also observed 
at 1.6 GHz. The catalogues, radio images and full de- 
scription of the observations can be found in the follow- 
ing papers: the Pearson & Readhead VLBI Survey (PR; 
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[14 1 lithe first Caltech-Jodrell Bank VLBI Survey (CJ1; 
[ 1§ |l|, 0); the second Caltech-Jodrell Bank VLBI Sur- 
vey (CJ2; jl8[ |ll|). From these samples a large sub-set 
of 300 sources was selected on the basis of a "flat" in- 
tegrated radio spectrum between 1.4 GHz and 5 GHz. 
In this context "flat" implies a source whose spectral in- 
dex a > —0.5, where the flux density at a frequency v 
is proportional to v a . The sample thus selected contains 
nearly all the strongest flat-spectrum radio sources in the 
sky above declination 35°. 

Flat-spectrum radio sources are best-suited to a 
millilcns search because they tend to be dominated by a 
single compact (sub-mas) "core" , which originates close 
to the massive central black hole in the nucleus of the 
background galaxy. The presence of two compact compo- 
nents (putative images) then provides a simple diagnostic 
of gravitational lensing by a point-like mass and genuine 
images will have well understood, simply-related, prop- 
erties || . However great care has to be taken in choosing 
lens candidates since compact radio sources have a range 
of intrinsic structures which can confuse the selection pro- 
cess. Guided by a series of simulations we decided only 
to look for cases of multiple imaging of the compact core. 
It is much harder to be certain of recognising lensing ef- 
fects on lower-brightness emission regions. We adopted 
the following quantitative selection criteria based on the 
Gaussian models listed in the original VLBI survey pa- 
pers: 

A candidate has to have two or more compact (< 1 
mas) components and the secondary component should 
be smaller than the primary. Because of the surface- 
brightness conserving property of lensing, weaker images 
are smaller images. However, conservative allowance was 
made for Gaussian model-fitting uncertainties by accept- 
ing some candidates in which the secondary had an area 
up to twice, and in a few cases three times, that of the 
primary. The separation of the primary and secondary 
components, 1.5 < A6 < 50 mas. This is set by the 
VLBI resolution and field-of-view limitations and deter- 
mines the mass range to which the search is sensitive. 
The upper limit on the separation was set by the prac- 
ticalities of map-making at the time the VLBI surveys 
were made. The ratio of the flux densities of the core of 
the primary component and the secondary component is 
< 40:1. At this level secondaries are detected with high 
signal-to-noise ratio. 

Fifty lens candidates were selected of which six could 
immediately be ruled out from other published VLBI 
data. 



III. SCRUTINY OF LENS CANDIDATES 

Since the brightness contrast between the core and 
all other emission regions in a radio source increases 
with increasing frequency, the simplest way to determine 
whether two or more of the components are compact, 
and have properties appropriate for lensed images, is to 
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FIG. 1: Upper) VLBI map of a millilens candidate 0740+768 
at 5 GHz and 1 milliarcsec resolution; the separation be- 
tween the two components is 3.2 mas. Lower) VLBI map of 
0740+768 at 22 GHz and 0.3 milliarcsec resolution; the two 
components, barely resolved at 5 GHz, have markedly dif- 
ferent surface brightnesses. Rather than a characteristic lens- 
ing morphology 0740+768 displays the asymmetric "core-jet" 
structure commonly found in flat-spectrum radio sources. 



observe at a higher frequency and higher resolution than 
the original 5 GHz survey maps. The candidates were 
therefore all observed at 15 GHz using the Very Long 
Baseline Array (VLBA). A minority of the candidates 
were also observed at one or more of the frequencies 1.6, 
8.4 and 22 GHz. The observations and the radio maps 
will be described elsewhere. 

Many of the candidates turned out to be core-jet 
sources in which the jet is faint but contains a "hot- 
spot". The VLBA 15- and 22-GHz maps reveal the 
bright core and resolve the lower-brightness hot-spot. An 
example of such a candidate is shown in Figure |. The 
5 GHz map of 0740+768 shows a simple unresolved dou- 
ble source but the higher resolution 22-GHz map clearly 
shows that the compact core lies at the eastern (left) end 
with the other "jet-like" emission having markedly lower- 
brightness. Other candidates turned out to be Com- 
pact Symmetric Objects (CSOs) — an early (< 10 4 years) 
high-luminosity phase of large double radio sources . 
At high resolution CSOs often show a single central core 
straddled by a pair of outer "lobes" of significantly lower 
surface brightness; not only is there no detailed corre- 
spondence between sub-structures in the lobes but also 
their parities are incorrect for lensed images. 

Component spectra can also be used to rule out can- 
didates. For lensing on the present scales the time delay 
between the light paths ranges from a few seconds to a 
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FIG. 2: Limit on fico a t the 95% confidence level from the 
failure to detect gravitational lensing amongst 300 compact 
sources with image separations in the range 1.5-50 milliarc- 
sec. 



few hours, which is much shorter than the timescales be- 
tween sets of VLBI observations. Lensed images must, 
therefore, brighten and fade together. A strong corollary 
is that lensed images must have identical radio spectra 
or, equivalently, their flux ratio must not be frequency- 
dependent. This argument should be used with care if 
one or both of the components is well resolved but we 
ruled out candidates with compact components if the 
flux ratios at different frequencies were inconsistent by 
>30%. 

Component motion is also an excellent discriminant. 
If the lens moves across the line of sight the relative sep- 
aration of the images will change. For example a lens at 
z ~ 0.5, moving transversely with v ~ 1000 km s _1 , will 
produce a relative image motion ~ 1 x 10 -5 mas yr _1 . 
This is much too small to be measurable with VLBI and 
thus, if relative motion is detected, the secondary cannot 
be a lensed image. In a separate follow-up programme 
for the PR and C J VLBI surveys most of the sources have 
been observed more than twice over a period of several 
years Currently over 30 of our 50 candidates are 

known to show component motions greater than ~ 0.1 
mas yr _1 and hence can be ruled out by this means as 
well. 

All 50 candidates were ruled out, often for a combina- 
tion of the reasons cited above, and we now use this null 
result to set quantitative limits on supermassive f2co for 
objects uniformly distributed throughout the universe. 



IV. LIMITS ON fico 



We used the detection volume method introduced by 
Nemiroff and developed in detail for the specific case 
of a VLBI search for point masses by KKB. We will, 
therefore, only make a few points specific to our search. 
The most important observational limit is the flux ra- 
tio R(< 40 : 1) which sets the maximum source-lens 
impact parameter and hence defines the basic detection 
volume. The usual definition of "strong" lensing || refers 
to images with a flux ratio of < 7 : 1; our < 40 : 1 
flux ratio therefore corresponds to a configuration with 
a larger lens-to-source "impact parameter" and hence 
the cross-section for lensing is almost five times that nor- 
mally assumed for lensing calculations. Each source is 
observed with an approximately fixed angular resolution 
and hence minimum image separation S(= 1.5) mas, and 
limited field-of-view A(= 50) mas, both of which act to 
truncate the detection volume. In effect 8 and A define 
the lower and upper limits of a mass range for which 
the truncation is small. Our search is sensitive to the 
mass range ~ 10 6 to ~ 10 8 M Q . For observations of a 
number of sources the individual detection volumes for 
each source can be summed to give a total detection vol- 
ume. In the calculation of these volumes one needs to 
know the source redshifts and these are available for 270 
out of the 300 sources in our sample; the mean redshift 
is 1.30 and 20% of the sources have redshifts > 2. For 
the remaining 30 sources we have assumed a distribution 
in redshift similar to that exhibited by the other 270. 
For the calculation of the angular size distances we took 
Hq = 65 km s _1 Mpc -1 and, for direct comparison with 
KKB's result we assumed an Einstein-de Sitter Universe 
(ftM = l;O A = 0). 

Our null result allows an upper limit to be placed on 
f^co an d Figure || shows the limit at the 95% confidence 
level as a function of CO mass. The limit is approxi- 
mately constant from ~ 10 6 to ~ 10 s M Q and in this 
mass range fico < 0.013. If Oco was equal to 0.013 the 
"expected" number of lenses in the total detection vol- 
ume associated with our 300 sources would be three and 
the fact that none are detected allows us to reject this hy- 
pothesis with 95% confidence. If f^co was equal to 0.026 
then ~ 6 lenses would be expected and we can reject this 
higher mass density with 99.7% confidence. 

These limits are conservative because gravitational 
lensing increases the observed flux density of a back- 
ground source and hence lensed sources are drawn from 
a fainter source population than the unlensed sources; 
a flux-limited survey will contain more lenses than ex- 
pected. Our null result therefore corresponds formally to 
a stronger limit on ilco but the "magnification bias" as- 
sociated with flat spectrum radio sources, especially for 
lens systems with high flux ratios, is of order unity p3] 
and we have therefore ignored this effect. 
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V. DISCUSSION 

KKB derived a limit on n co (10 7 -10 9 M ) < 0.4 
(99.7% confidence); this corresponds to f2 co < 0.2 (95% 
confidence). The lower mass limit is larger than ours 
since they took S ~ 4 mas compared with our 1.5 mas. 
The upper mass limit derived by KKB is high since it 
is based on A = 500 mas which is, in our view, opti- 
mistic; we used A = 50 mas. Our null result implies 
^co J; 0.013 (95% confidence) in the mass range ~ 10 6 
to ~ 10 8 Mo which is about 15 times more stringent than 
KKB's. The improvement arises from a combination of 
more sources (300 cf. 48), larger R (< 40 : 1 cf. < 20 : 1) 
and hence larger lensing cross-section and a higher mean 
redshift for the parent sample (1.30 cf. 0.89). Our result 
allows the following conclusions to be drawn. 

Uniformly distributed CO in the mass range ~ 10 6 
to ~ 10 s M & do not make up more than ~ 1% of the 
closure density fltotai — 1 which is strongly supported 
by the latest measurements of the angular spectrum of 
the Cosmic Microwave Background p3| . Similarly such 
CO do not make up more than ~ 3% of the Dark Matter 
density £Idm ~ 0.3 favoured by current observations. 
The favoured value of the baryon density from Big-Bang 
Nucleosynthesis is n b h 2 = 0.019 ± 0.002 @. Taking a 
plausible value for h (0.65) implies fl b = 0.045 ± 0.005 
thus our 95% confidence limit implies that uniformly- 
distributed Jeans-mass CO do not make up more than 



about one third of fib- 

A large population of uniformly-distributed ~ 10 6 ' 5 
Mq black holes is required by the cosmogonic model de- 
veloped by Gnedin & Ostriker S and Gnedin, Ostriker 
& Rees H . This model predicts that up to 5% of high 
redshift sources should be detectably lensed by such a 
population i.e. there should be up to 15 lenses in our 
sample. They are not observed and hence the model can 
be ruled out. 

Perhaps the next interesting limit is flco < 0.005 
which would constrain the contribution of supermassive 
CO to be no more than the baryonic contribution of 
presently observable stars and galaxies. To reach this 
limit about 1000 sources, with the same redshift distri- 
bution as the present sample, would have to be studied; 
this would be a time-consuming, but relatively straight- 
forward, task. 

We thank Silke Britzen for providing results prior to 
publication. The VLBA is the Very Long Baseline Array 
and is operated by the National Radio Astronomy Obser- 
vatory which is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associated 
Universities, Inc. 

Note added in proof. When this work was almost com- 
pleted we became aware of the millilensing results of 
Nemiroff and collaborators (preceding PR Letter) using 
gamma ray bursts. Their work has produced limits on 
the CO abundance which are similar to ours. 
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